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Abstract: We propose a new holographic model of Josephson junctions (and networks 
thereof) based on designer multi-gravity, namely multi-(super)gravity theories on products 
of distinct asymptotically AdS spacetimes coupled by mixed boundary conditions. We 
present a simple model of a Josephson junction (JJ) that exhibits the well-known current- 
phase sine relation of JJs. In one- dimensional chains of holographic superconductors we 
find that the Cooper-pair condensates are described by a discretized Schrodinger-type equa- 
tion. Such non-integrable equations, which have been studied extensively in the past in 
condensed matter and optics applications, are known to exhibit complex behavior that in- 
cludes periodic and quasiperiodic solutions, chaotic dynamics, soliton and kink solutions. 
In our setup these solutions translate to holographic configurations of strongly-coupled 
superconductors in networks with weak site-to-site interactions that exhibit interesting 
patterns of modulated superconductivity. In a continuum limit our equations reduce to 
generalizations of the Gross-Pitaevskii equation. We comment on the many possible ex- 
tensions and applications of this new approach. 
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1. Introduction 



There has been recent interest in the potential applicability of AdS/CFT- inspired methods 
to traditional condensed matter problems, which are not amenable to a weak coupling 
quasiparticle description. This has led to the formulation and study of a large variety 
of models in classical gravity which aspire to capture holographically the characteristic 
features of some condensed matter system. A notable example is [|l], Q which provides 
a simple gravity dual for an s-wave superconductor. More recent developments in this 
subject are reviewed in p, |, [|. Ultimately one hopes that holographic techniques will 
provide a new efficient description of high-T c superconductors that goes well beyond the 
BCS theory. 

Many technological applications of superconductors and superconducting devices in- 
volve Josephson junctions (JJs). The basic junction consists of two superconductors sep- 
arated by a weak link. The precise type of JJ depends on the specifics of the constituent 
superconductors and the nature of the link. The link can be an insulator (SIS junctions), 
a normal conductor (SNS junctions) or another superconductor. The coupled supercon- 
ductors can be of the same or different type. For example, one can consider sis, sld, or 
did junctions (s denoting an s-wave superconductor, d a d-wave superconductor and / an 
insulator). The properties of these junctions can be considerably different. For instance, 
quantum tunneling in conventional SNS and SIS junctions implies a current I across the 
link, even in the absence of external voltage, which depends on the phase difference $ of 
the condensates of the two superconductors in the following way |J 

I = /max Shit? . (1.1) 

This simple sinusoidal relation can be substantially different in other (more general) types 
of junctions (see for example ]7|, ^|, |S], |TU|, |TTJ). 



Another reason to be interested in JJs is the nature of high-T c superconductivity itself. 
Many high-T c superconductors enjoy a layered structure |12|] that can be viewed as a natural 



stack of atomic scale intrinsic JJs 1T3J with interlayer spacing of about 15.5 A. In fact, 



such high quality SIS-type intrinsic JJs can be fabricated [[14], [15] and pose as attractive 
candidates of cryoelectronics. 

Therefore, assuming holography can provide a new window to high-T c superconduc- 
tor physics, there is an obvious interest to construct the holographic dual of Josephson 
junctions and more generally the dual of Josephson junction networks (JJNs). 

A first step towards the construction of a holographic SNS junction has been taken 
recently in Ref. fl6j (see also |L7], [18|, |l^, for related setups). In this approach 
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one is looking for solutions of the equations of motion of the standard holographic su- 
perconductor setups that are inhomogeneous across one of the field theory directions (the 
direction along which the SNS stack is arranged). In the present paper we will propose 
a distinctly different way to construct a holographic Josephson junction (one that is not 
necessarily restricted to SNS types). 

The basic idea is to view each of the superconductors that compose the junction as 
a separate supergravity (or superstring) theory on its own asymptotically AdS spacetime 
and to model the weak link between them as a mixed boundary condition that relates the 
boundary conditions of the condensing symmetry-breaking field on one spacetime to the 
boundary conditions of the condensing symmetry-breaking field on the other spacetime. 
This may look like a contrived operation on the gravitational side but it is a rather natural 
one from the perspective of a dual large- N quantum field theory. On the field theory side 
this operation amounts to a multi-trace deformation that involves products of single-trace 
operators from both theories. This multi-trace interaction is the only term that mediates 
interactions between the two theories. The single-trace operators that appear in these 
interactions are charged under the broken U(l) symmetries and thus mimic naturally the 
charge quantum tunneling effects that are present in a JJ. 

Using multi-trace interactions to model the tunneling effects is natural in a regime 
where the mass scales of the modes that mediate the interactions between the two super- 
conductor field theories in the full system are large compared to the typical energy scales 
that we consider. In that case one can integrate out these higher mass modes to obtain an 
effective theory at low energies. The separate gauge invariance in each of the two boundary 
theories implies that the effective interactions between them can only be of the multi-trace 
type. Abstracting from this picture the main features of the bulk description we will pro- 
ceed to employ them freely in more generic situations where explicit knowledge about the 
boundary description is very limited or altogether absent. 

The precise ingredients of our construction are presented in section 0. A simple char- 
acteristic example of a holographic JJ at zero temperature is discussed in section [5| We 
show that the standard sine expression for the Josephson current ( |1 . 1|) is naturally repro- 
duced in this model in a few lines and determine J max explicitly in terms of the parameters 
of the system. We briefly comment on the extensions and modifications that can alter 
this standard current-phase relation. We also discuss how this system differs from a typi- 
cal Josephson junction and what kind of extensions can be used to describe more typical 
Josephson junctions. 

Another appealing feature of the above approach is its versatility in describing very 
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general configurations of Josephson junction networks. JJNs is a much studied subject 
with diverse applications. The great wealth of possibilities that they pose and the great 
reliability of the fabrication technologies developed for their construction makes them a 
prototype of complex physical systems that exhibits a variety of interesting physical be- 
haviors. Among their many applications: they are used widely as microwave sources (see 
e.g. ||23||), they provide controllable settings to investigate properties of granular or high- 
T c superconductors |p4]| , they are frequently used as model analogs of physical systems 
with complex dynamics. For instance, they have been used to model biologically realistic 
neurons [ f25|j . 

As a first simple application of our proposal in this direction we consider in section 
|] a holographic Josephson junction array that can also be viewed as a special case of the 
honey-comb network of Ref. p6fl . Using the gravitational description we find that the 
Cooper-pair condensates are described by a discretized Schrodinger-type equation, which 
has been studied extensively in the past (for a special value of one of our parameters) in 
radically different condensed matter and optics applications p7| . Using well known facts 
about this equation we show that the system in question exhibits complex behavior that in- 
cludes periodic and quasiperiodic solutions, chaotic dynamics, solitons and kink solutions. 
In our setup these solutions translate to one-dimensional configurations of holographic su- 
perconductor layers that exhibit interesting patterns of modulated superconductivity. In 
a continuum limit, the discretized Schrodinger equation becomes naturally a generaliza- 
tion of the Gross-Pitaevskii (GP) equation, a well-known long-wavelength description of 
superfluids. In this limit we recover some of the previously discrete solutions analytically. 

We conclude in section |5| with an outline of further possible extensions and applications. 



2. Networks of large- N QFTs and designer multigravity 

2.1 Networks of large-iV QFTs 

Consider a set (network) of k <i-dimensional QFTs with a large-iV limit — for example, 
k potentially different large- N conformal field theories (CFTs). We will label each CFT 
by an index i (i = 1, . . . , k). Equivalently, i is an index that labels a site (vertex) in our 
network. The links of the network are provided by interactions coupling the CFTj's with 
each other. The only kind of coupling respecting the individual gauge structure of each 
site is one mediated by multi-trace operators, so this is the only kind that we will consider 
here. For example, if Oi is a single-trace operator in the CFTj, then a double-trace link 
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vertex 




(a) 



o 




Figure 1: Examples of simple networks. In our context, these graphs will represent large-N 
CFTs coupled by double-trace interactions. 

between the CFTj and the CFTj corresponds to a Lagrangian interaction of the form 



Assuming that the large- N scaling of the single-trace operators Oi goes like O(N), the 
double-trace interaction (|2.1|) respects the large- N expansion when h is taken to scale as 
a constant, namely h ~ O(N ). Similarly, an £-trace coupling that involves the product 



It will be useful to set up a convenient notation to denote graphically networks con- 
structed in this way. For CFTs coupled by double-trace interactions, like that in eq. (|2.1| ), 
we will denote the corresponding link by a single line. Such a link is undirected and joins 
two different CFTs, or circles back to the same CFT. The latter denotes that the corre- 
sponding CFT has itself a double-trace interaction turned on. Hence, Fig. |I](b) exhibits a 
simple network of two CFTs, call them CFTi and CFT2, with the Lagrangian interaction 



Similarly, Fig. |I|(c) exhibits nine CFTs coupled pairwise by double-trace interactions in all 



6C = h OiOj . 



(2.1) 



flf=i ®is should have a coefficient that scales as 0(N 2 e ). 



5C = hnOf + h 12 O x 2 + h 22 Ol . 



(2.2) 



possible ways. This exampl< 
using the replica trick p9| . 



.e, analyzed in |28[, plays a role in quenched disorder calculations 
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(c) 



Figure 2: Figures (a) and (b) exhibit a triple-trace and four-trace link respectively. In order to 
make the notation more transparent we have denoted a multi-trace link using a double-line notation 
reserving the single-line notation for the simpler double-trace link. Accordingly, Figure (c) exhibits 
four CFTs linked pairwise with two double-trace interactions whereas Figure (b) exhibits a four- 
trace link. 



Networks formulated in this way can have different types of links. We may use a 
double-line notation to denote the more general possibility of links mediated by multi-trace 
interactions. Figs. |](a) and ||(b) exhibit a triple-trace and four-trace coupling respectively. 
In network literature the corresponding graphs are sometimes called hypergraphs. The 
use of such more general couplings opens up many interesting possibilities. Most of our 
discussion in this paper will be focused, however, on networks with double-trace couplings 
only. 

Being comprised by sites that correspond to interacting QFTs the above networks have 
in general a rich and complicated internal structure. This structure can evolve in time, 
vary in space, or vary from site to site. Moreover, under renormalization group flow the 
network graphs may change with the appearance of new links or even new sites. Indeed, it 
is well-known that multi-trace couplings are naturally generated under renormalization (see 
for instance PS|). Hence, renormalization effects in field theory can affect the quantitative 
features of the links. They can also change the number of sites in the following way. Many 
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supersymmetric QFTs, like the M = 4 super- Yang-Mills theory, have a Coulomb branch. 
At generic points of this branch the original gauge group is Higgsed into product gauge 
groups and in the far infrared one is left with a product of QFTs giving rise to a network 
with more sites than those of the UV theory. This is an example of an RG flow with 
a different number of sites in the UV and IR. In what follows, when we draw a graph 
representing a network we will implicitly assume that this description refers to the bare 
UV Lagrangian of the corresponding theory. 



2.2 Networks of asymptotically AdS spacetimes 

There are situations where a large- iV QFT (typically at strong coupling) has a dual descrip- 
tion in terms of a supergravity theory on an asymptotically AdS background. Accordingly, 
a network of such QFTs has a dual description as a multi-gravity network where each site 
is some supergravity theory on an asymptotically AdS background and each link is a mixed 



boundary condition for supergravity fields residing on different space-times [|30], |31], |28| (see 



p2| for a stringy setup that involves multi-string theory networks and |33[ for a brief review 
of the main idea and its implications for massive multi-gravity). Let us recall how such a 
multi-gravity network comes about in the AdS/CFT correspondence. For clarity, we will 
restrict to the case where the boundary QFTs are conformal. 

Before the addition of multi-trace links, the CFT^s are independent field theories that 
do not talk to each other and the dual gravitational theory is a direct product of super- 
gravity theories on product AdS space-times of the form n?-=i t^*^ ® -M«]> where Aii is 
the internal manifold of the spacetime of the dual of CFTj. 

Assume that the single-trace operators Oj stated above are scalar operators 1 with 
scaling dimension Aj. The AdS/CFT correspondence maps each Oi to a dual scalar field 
ifi with AdS asymptotics 

<Pi*% + . + + • (2.3) 

i i 

r,i is the radial distance in the 2-th AdS space that corresponds to CFTj. We use conventions 
where the i-th AdS boundary lies at r« — > oo. In later applications, we will assume Aj < |, 
in which case the a^s should be interpreted as the vacuum expectation values (VEVs) of 
the dual operators and the /%'s as the sources. 



-"-It is not necessary to restrict ourselves to scalar operators. We will make this assumption here for 
reasons of simplicity and concreteness. In fact, some of the applications of this framework that we will 
propose later also involve vector operators. 
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In this framework, and to leading order in the 1/N expansion, it is well known that a 
multi-trace coupling in the boundary CFT 



5C = W(0 1 ,...,O k ) (2.4) 

imposes mixed boundary conditions to the asymptotic coefficients cij, of the form 

P l = d l W(a 1 ,...,a k ) . (2.5) 

Combined with requirements of regularity these boundary conditions fix completely the 
profile of the bulk solution. The inter-theory coupling induced by the relations (|2.5|) leads 
to a network of supergravity theories that can exhibit interesting collective phenomena. 
We will explore these phenomena in section |j. 

2.3 Designer multigravity 

Let us recall how multi-trace deformations affect the profile of the bulk asymptotically 
AdS solution in a single theory; a subject that usually goes under the title of designer 
gravity ||34|| . In this paper we want to consider the natural extension of this framework to 



multi-AdS spaces, which we will suggestively call 'designer multigravity'. 
As one of the simplest illustrations of the idea consider the 'network' 




(2.6) 



which contains a single vertex with a self-adjoining link. The graph Q2.6|) exhibits a double- 
trace link, but we can equally well consider any multi-trace link. This theory has been 
the main focus of most previous investigations of multi-trace interactions in the AdS / CFT 
correspondence and designer gravity. For self-completeness and in order to set the notation, 
we briefly review some of the most pertinent properties of this system. 

At the single site of this network resides a <i-dimensional large- N CFT with a super- 
gravity dual. Assume that the CFT has a single-trace complex scalar operator O with 
scaling dimension A. This operator is dual to a charged bulk scalar field (p. To further 
simplify the discussion we will also assume that we can consistently reduce the dynamics 
of the dual bulk supergravity to a (d + l)-dimensional Einstein- Abelian Higgs model of the 
form 

1 



Sbuik = I d d+l x 



R - -G(\^\)F 2 - (VM) 2 - J(M) (W - qAf - V(M) 



(2.7) 
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where 9 is the phase of if, namely (p = Ifle 10 , and G, J, V are model-dependent functions 
of \ip\. A is an Abelian gauge field in the bulk and F its field strength, q is the U(l) charge 
of the field (p. 

We are interested in asymptotically AdS^+i solutions of this system. In units where 
the AdS radius is set to one, the potential V(ip) has the small-y? expansion 

V(\<p\) = -d(d - 1) + m 2 M 2 + 0(M 4 ) + . . . . (2.8) 

When the mass m lies within the range 

d 2 

m 2 BF < m 2 < m 2 BF + 1 , m% F — — — (2.9) 
the dual operator O can have two possible scaling dimensions 



d d 2 

A ± = 2±VT + m2 - (2 - 10) 



We will assume that our theory lies in the window fl2.9|) and will pick O to have the smaller 
scaling dimension A = A_ that satisfies the inequality < A_ < | (the lower bound in 
this inequality is the standard unitarity bound in field theory). 

Near the asymptotic boundary, r — > oo, the metric is that of AdS^+i 

dv 2 

ds 2 -^dx^dxt* + — (2.11) 



r 

and the scalar field <p exhibits two independent branches 

<P*£ + - + -L + - ■ ( 2 - 12 ) 

The boundary condition 

>-%• 

for a generic smooth function W(a) corresponds at the boundary CFT to the multi-trace 



deformation 35 



5C = W(0) . (2.14) 

When we solve the bulk equations of motion we are looking for solutions that respect 
the boundary conditions ( |2.13|) . In addition, we require that these solutions are every- 



where regular. It turns out that regularity imposes an extra constraint on the asymptotic 
coefficients a, (3. We will denote this additional relation as 
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where W is a function determined by the specific dynamics of the theory, e.g. the details 
of the scalar potential V(|<^|) in the bulk action ( p.7|) . 

Combining the boundary condition (|2.13|) and the regularity condition ( p. 15 ) we find 



that a and (3 are completely fixed and determined as a solution to the equation 

dV ~ 

— = 0, V(a) = W(a)+W(a) (2.16) 
da 

which therefore can be viewed as an extremum of the function V. 

For example, for boost invariant, planar configurations with vanishing gauge field 

dv 2 

ds 2 = r 2 i-dt 2 + dxidx 1 ) + — -r , <p = <p(r) , A = (2.17) 

gir) 

one finds a solution with an acceptable naked singularity at r = that has 

W(a) = s^\a\ d/A . (2.18) 
d 

The existence of this solution and the precise value of the parameter s depends on the 
details of the bulk potential V((p) (see |3(| for additional information). For most potentials, 
s turns out to be a positive number. The energy density of this solution is 

A = ( d _2A)V. (2.19) 

One can also consider analogous solutions with spherical topology. It is possible to 
show that the energy fl2.19|) , ( [2.18] ) provides a lower energy bound for all of these solutions 



[ 36| . Moreover, in the case of vanishing gauge field, one can argue [3£| that 



• the theory with boundary conditions f3 = ^ has a stable ground state provided the 
function V has a global minimum V m in, and 

• that the minimum energy solution is the spherical soliton associated with V m j n . 

Note that because of the presence of W, it is possible to have a stable ground state even 
for functions W that have no minimum. In general, this minimum involves a condensate 
of the charged scalar that higgses the corresponding U(l). 

What we have discussed so far applies to the case of zero temperature. By study- 
ing hairy black holes in the bulk it is possible to generalize the discussion to non-zero 



temperature |37| . It is also possible to consider the case of non- vanishing charge density 



[38 



It is not hard to generalize this discussion to the case of multiple CFTs coupled together 
by a multi-trace interaction of the general form W(0i, . . . , Ok). In this case (and to leading 
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order in the 1/N expansion), the bulk equations of motion are the same as before in each 
(super)gravity member of the product, but the boundary conditions change. For example, 
the bulk boost-invariant planar solutions in the i-th spacetime still retain the form ( |2.17 ) 
and each (3{ is still given by an equation of the form 

ft = (2.20) 

with 

W l (a i ) = s^\a l f^ . (2.21) 

The new ingredient, responsible for the coupling between different AdS theories, lies in the 
mixed boundary conditions ( |2.5| ). The analog of equation ( p. 16 ) that determines the VEVs 
ai is 

&V(a 1 ,a 2 ,... 1 a>l =Q (2 22) 



da 
with 



A: 



V = WK, ...,«*) + £ Wi(a<) • (2.23) 

i=i 

In general, these equations lead to a non-linear discrete system which can exhibit intricate 
behavior. This behavior includes solutions that can be periodic, quasi-periodic, chaotic 
or even soliton-like with energy pinned around a central site. Examples of each of these 
behaviors will be discussed in section f|. 

2.4 A holographic superconductor with vanishing charge density 

A novel type of holographic superconductor with symmetry breaking induced by double- 
trace deformations was recently proposed in ||38|| . Since this setup will provide the basic 
building block of the discussion to come, it will be beneficial to recall some of its main 
properties. 

Returning to the single-site example of ( |2.6D consider the case of a double-trace defor- 
mation 

W = g\0\ 2 (2.24) 

implemented with the use of a single-trace operator O with scaling dimension A. The 
function V in (p. 16 ) becomes in this case (for boost-invariant planar solutions) 

V{a)=g\a\ 2 + S -\a\\ 6=±. (2.25) 

The extrema of this function obey the algebraic equation 

«(^ + ||a|^ 2 ) = • (2.26) 
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Assuming s > 0, there is one or two possible solutions to this equation depending on 
the sign of the constant g. For g > 0, the only solution is a = 0, which is a minimum and 

does not exhibit any condensate of the field (p. For g < 0, there are two possible solutions 

i 

« = , a e = f-^A S " e w . (2.27) 

The first one, cto, is a local maximum of V and therefore an unstable vacuum of the system. 
The other solutions, aig, are degenerate stable minima labeled by an angular variable 9. 
The non- vanishing condensate of the charged scalar field ip in this case higgses the U(l) 
gauge symmetry and leads to a new type of holographic superconductor. 



3. A holographic model of Josephson junctions 
3.1 The setup 

A Josephson junction consists of two superconductors separated by a link that mediates 
weak interactions between them (several possibilities for the weak link were reviewed in 
the introduction). In analogy, consider a setup where two holographic superconductors 
(each described by a gravitational theory on an asymptotically AdS space) are interacting 
weakly via mixed boundary conditions on the boundary. 2 From the gauge theory point 
of view two initially separate large- N gauge theories are brought into contact via multi- 
trace interactions. As we discussed in the introduction, in certain cases one can think of the 
multi-trace interactions as an effective description below the typical mass gap associated to 
the actual interactions between the two systems. In a real Josephson junction these would 
be the weak tunneling interactions across the material in between the superconductors. 
As a concrete illustration of the idea consider a pair of two identical holographic super- 



conductors of the type described in the previous subsection |2. 4 The parameters g, s, 5, d, A, 
which will be treated here as phenomenological parameters of the model, are chosen to be 
common in the two systems. Each of them has a charged scalar field 0{ [i = 1, 2) and a 
corresponding dual complex scalar field tfi with asymptotics (|2.12|) . We can think of Oi 
as the 'Cooper pair' operator in the holographic superconductor CFT,. Holographically, 
the VEV of these operators are given by the leading branch coefficients aij in the asymp- 
totic expansion (|2.12|). We assume that both operators have the same scaling dimension 



d 
2 

AdS spacetime. 



A < | and employ the alternative quantization for the dual scalar fields (fi in the bulk 



2 We will soon discuss how similar this system is to the typical Josephson junction considered in the 
laboratory and the extensions it suggests. 
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In order to mediate an interaction that exchanges charge between the two theories we 
couple them via a double-trace interaction of the form 



W(O x , 2 ) = h (e i ^0 1 0\ + e-^0\0 2 ] . {3.1) 



h is a real number controlling the strength of the interaction and $ is an angular variable. 
The latter is an additional tunable parameter of the interaction whose physical meaning 
will become clear in a moment. 

In the bulk this coupling implies mixed boundary conditions. One can determine 
the VEVs of the dual operators Oi by solving the scalar-gravity equations of motion. 
Equivalently, one can extremize the function V (see eqs. ( |2. 22| ) , (|2.23|) ). In the case at 
hand 

2 

V{a u a 2 ) = (g\®i\ 2 + ^N 5 ) + h (^V^ + e~ i{> a* l a 2 ) . (3.2) 

i=l 

Extremizing with respect to a\ and a 2 we obtain two algebraic equations 

s 

ga\ + he~ l ^a 2 + -ai\a\\ s ~ 2 = , 

g^ + he^^ + ^a^] 5 - 2 = (3.3) 

which determine the a-j's uniquely in terms of the parameters g, h, s and 5 up to an overall 
phase common in a± and a 2 . The relative phase between a\ and a 2 is fixed in terms of $ 

tf 12 = $ 2 - #i = # mod 7t . (3.4) 

For instance, when 5 = 4 3 the solutions are 

(1) ai = 0, (2) ± \ ai \ 2 = -(±h-g) , (3)± \a 1 \ 2 = --(g±y/g 2 -4hA . (3.5) 

s s \ J 

In all cases 

a 2 = -/i"V 1, ^ + ||ai| 2 )ai. (3.6) 

Which solutions survive in cases (2) and (3) depends on the specific range of the parameters 
s,g,h. Assuming s > we obtain 

• \h\ < g: vacuum (1), 

• -\h\ < g < \h\: vacua (1), (2) sgn(/l ), 



3 This is consistent with Ae (f — 1? |) iff A < 1, which is possible only if d — 2,3 (the main cases of 
interest in this paper). 
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• —2\h\ < g < —\h\: vacua (1), (2)±, 

• g < —2\h\: vacua (1), (2)±, (3)±. 

Diagonalizing the Hessian of the energy functional V (see eq. ( [2.19] )) we find that (1) 
has the eigenvalues 

2(p + h) , 2(g - h) . (3.7) 

Hence, when phases with a non-zero condensate exist (g < \h\) the vacuum (1) of the 
normal phase is unstable (one of the eigenvalues of the Hessian is negative). In that 
case there is always at least one superconducting vacuum that is stable. Such a phase 
provides a holographic model of two weakly interacting superconducting materials at zero 
temperature. 

As a simple illustration, by setting h = 1 we find that 

• when — 1 < g < 1, the only stable vacuum is (2) + , 

• when —2 < g < — 1, there are two stable degenerate vacua (2)±, and 

• when g < —2, the vacua (3)± are also unstable and the only stable vacua are again 
(2)±- 

3.2 Josephson current 

A characteristic feature of conventional Josephson junctions is a transverse supercurrent I 
which is related to the condensate phase difference $12 in the following way 

I = /max sin ^12 . (3.8) 

In this section we have considered a system of two sites that represents two infinitely 
thin layers of superconducting material at zero charge density coupled through a weak link 



expressed by an interaction of the form (|3.1|), which can be suggestively rewritten as 



W(Oi, 2 ) = h {fO x O\ + e-^0\0 2 ) =W E + W Jaxt (3.9) 
with the definition 

W E = hco8#(0 1 0l + 0l0 2 } , Wj emt = ihan-&(OiOl-Ol0 2 > \ . (3.10) 

We is an interaction that mediates no interlayer charge transfer. In contrast, yVj ext is an 
interaction based on the charge-transferring operator 

J esBt = i(0 1 0t-0\0 2 ) . (3.11) 



- 14 - 



The coupling A of this operator in Wj ea . t , namely 



W Jext = AJ ext , A = hsm$ (3.12) 

can be interpreted as a new interlayer background gauge potential component. 4 In this 
sense, for non- vanishing h, $ our two-site system lies in an external transverse gauge field 
and J ext is an externally imposed current. 

The total current running across the sites of the junction can be determined in the 
standard way from an infinitesimal relative U(l) gauge transformation of the action. For 
an interaction of the form ( ^.9[ ) the infinitesimal transformation 

8O1 = ieOx , 50 2 = -ie(D 2 (3.13) 

gives the Lagrangian variation 

5C = 2ihe (e^O.Ol - e~^0\0^ = e(J slte x - J slte 2 ) = 2eJ tot . (3.14) 

The second equality is the discretized version of the gradient dJ across the interlayer 
direction. In the third equality we used charge conservation to set J site % = — J site 2 = Jtot- 
Consequently, 

J tot = h (fO x O\ - e- w O\0^j . (3.15) 

Hence, in the vacuum governed by the algebraic equations ( pj.3[ ) one finds (at leading order 
in the 1/N expansion) 

J tot = h {e^onal - e- w ala 2 ) = . (3.16) 

This is precisely what one expects. Since our system is kept at zero charge density, in the 
equilibrium state charge cannot flow across the junction between the two sites. As a result, 
irrespective of the initial configuration, once the interaction ( |3.9| ) is turned on the system 
backreacts and evolves to a new vacuum, which can be conveniently determined with the 



holographic techniques of subsection |3.1| . In the new vacuum, which is characterized by 
the solutions of the algebraic equations ( |3.3j ), the condensate phase difference is $12 = — $ 
and the magnitude of the condensate has adapted accordingly and in direct relation to the 
strength of the interlayer couplings h, 



4 This component of the gauge field is not an inherent quantity of the (d + l)-dimensional theories 
living in each site of our network. In general, a linear array of sites 'deconstructs' an extra spacetime 
dimension and quantities, like A, associated with this extra direction arise in the field theory space of the 
lower-dimensional multi-AdS/CFT network (quiver) as new interlayer interactions. 
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It is interesting to consider the VEV of the externally forced current operator J ext in 
the new vacuum. At leading order in the 1/N expansion 



(Jext) = i(ot\a 2 - aia* 2 ) . (3.17) 

Using eq. ( |3.6| ) we obtain 

(J ext ) = /max sin ■& = -J max sin # 12 (3.18) 

with 

7___ = 2h- 1 \ry,\ 2 (a 4- 



2/i- 1 |«i| 2 (<7 + Jail 2 ) • (3.19) 



In the algebraically simple case of 5 = 4 



s 

on the vacua (2)± and 



J max = ±- (±/i - (?) (3.20) 



Ah 



1 max 

S 



(3.21) 



on the vacua (3)± (interestingly in this case J max is independent of g). 

The vanishing of the total interlayer current J to t implies that the backreaction has 
created an equal and opposing 'Josephson current' Jj OS e P hson (due to the condensate phase 
difference), which cancels the contribution of the externally imposed (J ex t)- Specifically, 

J josephson {Jext) -^max sin (3.22) 



in agreement with the expected sine law (|3.8|). 

It is clear that the system we have just described is a peculiar Josephson junction 
unlike the typical Josephson junction commonly discussed in the literature and engineered 
in the laboratory. In contrast to our system in a typical junction the superconductor com- 
ponents have finite spatial thickness in the transverse junction direction, charge can flow 
in this direction, the condensate phase difference is dialed by choice (and not determined 
dynamically) and a Josephson current arises without having to apply a gauge field across 
the weak link. In section f| we will describe how such more conventional junctions can be 
engineered in our framework as simple extensions of the setup we have described in this 
subsection. 

3.3 Possible extensions and other current-phase relations 

We briefly comment on a few possible extensions of this framework that may ultimately 
lead to different current-phase relations. 
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A possible alternative is to consider the same superconductor models as above, but 
modify the double-trace interactions. For example, one could consider other types of 
double-trace interactions, e.g. instead of (or in addition to) ( |3.1|) a double-trace deformation 
of the form 

h{e w O x 2 + e-™0\0\) . (3.23) 



This will modify the equations (|3.3|) and the Josephson current (|3.22|) in an obvious fashion 



Another possibility is to consider higher multi-trace deformations. An example of a 
triple-trace interaction is 

h[e w O x O^ + e- iu, 0\0 2 2 ) . (3.24) 

Finally, it is interesting to consider coupling different types of asymptotically AdS 
theories. For example, one can try to couple the holographic s-wave superconductor of Ref. 
H to the holographic p-wave superconductor of Ref. [^]. In the s-wave superconductor it 
is a bulk complex scalar field, dual to a complex scalar operator, that condenses. In the 
p-wave superconductor there is an SU(2) gauge field in the bulk and the U{\) C SU(2) 
is broken by the condensation of the remaining two components of the SU (2) gauge field, 
which are dual to the corresponding components of an SU(2) current in the boundary 
theory. One can consider the possibility of scalar-current double-trace deformations on the 
boundary which presumably translate to mixed boundary conditions for the dual scalar 
and gauge field components. 



4. Holographic Josephson junction arrays 

A lot of theoretical and experimental work has been performed on Josephson junction 
arrays (and more generally networks of Josephson junctions). Some of the main motivations 
and results in this field were summarized in the introduction. A natural generalization of 
the holographic construction of the previous section can be used to model networks of 
holographic superconductors. 

In what follows we will concentrate on a simple network with the topology of a chain 
as depicted in Fig. ^. In this network M sites (labeled by an index n, each of them 



representing a d- dimensional holographic superconductor of the type of subsection pT4|) are 



linked by the mixed boundary conditions corresponding to the double-trace interactions 
W({0„}) = 9 °A + h (j»O n Ot +1 + e^OlO n+1 ) . (4.1) 
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Figure 3: A chain of holographic superconductors labeled by an index n £ Z. .Eac/i Zinfc denotes 
an interaction mediated on the field theory side by a double-trace deformation. 

For simplicity, we are assuming the same constants g, h, S for all sites and links. Ac- 
cordingly, the potential function V that we have to extremize in gravity to obtain the 
condensates a n is 

V({a n }) = J2(9\ a n\ 2 + ^\a n \ S + h(e^a n a* n+1 + e- i& a* n a n+1 )^ . (4.2) 

n 

The extrema of this function are sequences of complex numbers obeying recursion relations 
with rich features. The solutions can be periodic, quasiperiodic (chaotic) or solitonic, and 
provide interesting new examples of spatially modulated, namely lattice site-dependent, 
superconductivity. 

We will organize the discussion according to the number of boundary conditions. 
4.1 No boundaries 

In the case of no boundaries, the extremization of the potential function V (|4.2|) with 
respect to all the otj's gives the recursion relations 

ga n + h(e i ^a n - 1 + e-^a n+1 ) + -a n \a n \ 6 - 2 = , neZ. (4.3) 

Setting 

a n = e <rt V, ~9=l , ~s = { (4-4) 
we can recast ( |4.3|) into the form 

9<Pn + fn-\ + ¥n+l + T^nl^nl 5 " 2 = . (4.5) 

The generic solution of these equations is parameterized by two complex numbers 
(these could be for instance the values of <po, <pi at the vertices and 1). There are 
various ways to analyze this system of equations. In fact, such equations have appeared 
in the past in a variety of applications and discussions of discrete dynamical systems. A 
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notable example, that arises for the special case of 5 = 4, is that of the discrete non-linear 
Schrodinger (DNLS) equation 

-^VVi + 1>n-l + 1pn+l + ^n\^Pn\ 2 = . (4.6) 

This equation is a discretized version of the Schrodinger equation with a non-linear quartic 
potential. Setting 

V> n = ^ n e lEt (4.7) 

we recover our set of equations (|4.5| ) with E = g and 5 = 4. 

The DNLS equation has a long history (for an extensive review and references we refer 
the reader to |]27||). In solid state physics it appeared first in the context of the Holstein 



polaron model for molecular crystals |40| . In an optics context, DNLS describes wave 
motion in coupled nonlinear waveguides. In the basic nonlinear coupler model introduced 
in |41] two waveguides made of similar optical material are embedded in a different host 
material. DNLS, and generalizations, have also appeared in studies of nonlinear electrical 
lattices fl4"2f . 



As a concrete algebraically simple example, in the rest of this section we will concen- 
trate on the case of 5 = 4. Qualitatively similar results are expected for generic 5. The 
linear stability analysis is treated for generic 5 in appendix [A]. For most purposes we follow 
closely the analysis of ]27| which we recommend for additional details. 



4.1.1 The general structure of solutions 

The recursion relation ( [4.5|) (with the ansatz 5 = 4) 



gtp n + if n ^i + tp n+1 + -if n \if n \ 2 = (4.8) 
can be viewed as a four- dimensional mapping from C 2 — > C 2 . Using polar coordinates 

ip n = r n e ie - (4.9) 
we obtain the following two sets of equations 

r n+ i cos(A0 n+ i) + r n _i cos(A0 n ) = - (g + '^r^j r n , (4.10) 

r n+1 sin(A0 n+ i) - r n _i sin(A6» n _i) = , (4.11) 

where 

Ae n ee e n - n _! . (4.i2) 
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Equation ( |4.11D is equivalent to the conservation of current 



J = TnTn-i sin(A0 n ) . (4.13) 
It is convenient to introduce the real-valued variables 

x n == (f* n (f n -i + PnVn-i = 2r„r n _i cos(A6> n ) , (4.14a) 

Vn = i (VnVn-X ~ <Prti>n-\) = 2J ' (4.14b) 

Zn = \¥n\ 2 - \¥n-l\ 2 = T* ~ r\_ x . (4.14 C ) 

In terms of these variables the recursion equations (|4.8|) become 



( //- | (?/;„- z n ) ) / :„ ) . ( 1. 15a) 

z n+l + z n = l - Xn+ \ Xn , (4.15b) 
2 w n + z n 



w n = y/x* + z 2 n + 4 J 2 (4.15c) 

thus reducing our 4D map to a 2D map .M : IR 2 — ?• M 2 . 

This map depends on two parameters: (g, s). The dependence on J can be scaled away 
by setting 

x n — > 2Jx n , z n — > 2Jz n , 5 — > 2Js (4-16) 



so that w n = a/1 + x 2 n + z\. A linear stability analysis shows (see appendix |A] for details) 
that, depending on the precise parameters, there are both bounded and diverging solutions. 
In certain regimes, e.g. when 

3^9 2 (^ + 2 ) , \ m .f-a . 2(2 -g) 
s > , # < 2 , < \ip n \ < — — — (4.17) 

s s(o — 1) 



the solutions are regular and bounded and the Lyapunov exponent vanishes \ 2q . Recall 
that the original parameter s that appears in the AdS / CFT context ( |2.18| ) is positive, but 
s — I can be both positive or negative. Moreover, one can tune s freely by adding on the 
field theory dual quartic-trace interactions of the form \0\ 4 . 

The general structure of the space of solutions is organized by a hierarchy of periodic 



orbits surrounded by quasi-periodic orbits |27[]. The periodic orbits can be traced on the 



intersection of any two symmetry lines Sq = M u Sq } «Sf = M n Si, n = 0, 1, . . . ^4 
where the fundamental symmetry lines are defined as 

S : z = , (4.18a) 
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Si: x = -- ig+-(w + z)j {w + z) . (4.18b) 

For example, in the intersection of the lines So and Si one locates the fixed point (x = 
£*, z = 0) with 

x* = -\[9+ ^1 + ^) y/l + xl . (4.19) 

The linear stability of an orbit with period q is conveniently characterized by the value 
of Greene's residue [[45] 

1 



P 



2-Tr 



(4.20) 



f(n) 

l—S J v f 

where DAA is the linearization of the map Ai. The period orbit is linearly stable when 
< p < 1 (elliptic periodic orbit) and unstable when p > 1 (hyperbolic with reflection) or 
p < (hyperbolic). 

In the case of the fixed point (x*, 0) the residue is 



P 



+ | v 7 !^) (2 + 5 V / T+^!) • (4.21) 



For s = eq. (|4.19 ) has one root, x\ = 4^2, which is real when \g\ < 2. In that case, 
< p = 1 — jg 2 < 1, so one obtains an elliptic fixed point. This conclusion continues to 
hold for generic s > and \g\ < 2. 

In the regime of elliptic stability the fixed points (a;*, 0) form the largest basins of sta- 
bility among all elliptic orbits. These stable orbits, which include both periodic and quasi- 
periodic solutions, encircle the fixed point forming the main island on the map plane. The 
quasi-periodic orbits, which lie on closed curves (the Kolmogorov-Arnold-Moser (KAM) 
tori), densely fill the island. One can show that the map Ai is topologically equivalent 
to an area-preserving map ensuring the existence of such KAM-tori near the symmetric 



elliptic fixed points [ 46[ . 

An illustration of the main island of elliptic orbits around the fixed point (x*, 0) for 
g = 1.6, s = 0.1 can be found in plot (a) of Fig. |] Outside this island one finds regular 
quasiperiodic orbits and orbits that diverge. 

Interesting changes in the structure of the solution space can occur as we vary the 
parameters of the system, here (g,s). In particular, the residue of periodic orbits can 
change. When the residue changes from a positive value to a negative value then a tangent 
bifurcation occurs where an elliptic point converts into a hyperbolic point. Whenever the 
residue exceeds the value of one from below a stable elliptic orbit converts into an unstable 
hyperbolic point with reflection accompanied by the creation of two new stable elliptic 
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Figure 4: Plot (a) depicts in the {x,z)-plane the main island of elliptic orbits that develops 
around the elliptic fixed point (x*, 0) with ~ —1.59355 for g = 1.6, s = 0.1. Plot (b) depicts a 
characteristic example of period- doubling bifurcation for g = —0.62, s = 1. Orbits with different 
initial conditions are depicted with different colors. 



points. This kind of bifurcation is known as period-doubling bifurcation — a 1-period 
orbit (fixed point) converts into a 2-period orbit. The new elliptic orbits remain stable 
until another period-doubling bifurcation occurs. After a cascade of such bifurcations local 
chaos appears. 

In our system, the residue p in eq. (f4.21|) , is always less than one assuming s > 0, 



g > 0. In that case, only tangent bifurcations can occur and global chaos can arise through 
the so-called resonance overlap. Period-doubling bifurcation can instead occur when g < 0. 
Specifically, in the range g < — | < a new period-2 orbit is generated at the points 



x± = ±y4^-l, z = 0. (4.22) 

The map M. acts on these points by sending (x±,0) —> (x T ,0). The residue of the new 
solution is positive 

5 (>-?)■ < 4 - 23 > 

As we decrease g further the period-2 orbit loses its stability and a new period-doubling 
bifurcation occurs which gives rise to a period-4 orbit. This cascade terminates at a critical 



parameter (see 27] and references therein) 



- { s- \i -- iOoi j . <\ ~ i. -<>-.<. < i.2i) 

which is called the accumulation point. At that point local chaos appears. 
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A characteristic example of period-doubling bifurcation is depicted in plot (b) of Fig. 
|j. At g = —0.62, 5 = 1, the originally stable fixed point of plot (a) has become unstable 
and two new elliptic fixed points have been created giving rise to elliptic period-2 orbits. 

Besides the issue of linear stability, that was discussed above, one can also ask about 
the local and global thermodynamic stability of the above solutions. Local thermodynamic 
stability requires a positive definite Hessian of the multi-gravity energy functional ( [2.23D 
(see also (|2.19|) ). It would be interesting to examine the extent to which local thermody- 
namic stability is equivalent to linear stability. On the other hand, global thermodynamic 
stability implies that the more stable solutions have less energy. In the following subsec- 
tions |4.1.2| and [4.2| we will see that chains with a finite number of sites have a discrete set 
of solutions. In that case, the solutions with the minimum energy are thermodynamically 
favored. We hope to return to a more detailed examination of these issues in the future. 

4.1.2 Periodic boundary conditions 

From the above discussion it should be clear that, for given values of the parameters and 
period, one is left with at most a discrete finite set of solutions to the recursion relations 
( |4.(j| ) in the case of periodic boundary conditions. Indeed, in regimes that allow for elliptic 
periodic orbits the choice of a prescribed period picks the sequence of y? n 's in general 
uniquely. In other regimes of parameters periodic solutions do not even exist. Moreover, 
through period-doubling bifurcation it is interesting to note that it is possible to have 
periodic solutions where the cp^s arrange themselves in more than one different domains 
of values. 

4.2 One or two boundaries 

If there is a boundary, say at n = with n valued only on non-negative integers, the n = 



version of the equation (|4.<j|) is modified to 



gtpo + <pi+ 2<A)|<A)n =0 ■ (4.25) 

In that case, the whole solution is fixed by the choice of one parameter, for example (po. 
If the chain has finite size and there is also a second boundary, then the analog of ( [4.25 ) 



at the second boundary will fix <p as well and the solution will be discretely unique and 
expressed completely in terms of the parameters of the system g, s, h. A specific example 
of this situation appeared in the dimer case of section |3|. 

A special consequence of eq. (|4.25[ ) is the fact that <p n are all real-valued up to a 
common n-independent phase. Equivalently, the phases 9 n in (|4.9| ) are all equal modulo 
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7r. This property can be deduced by using the boundary equation (|4.25|) to compute the 
conserved (i.e. n- independent) current J ( }4.14b| ). One finds J = from which the above 
statement follows immediately. 

Taking care to satisfy the conditions from the boundary equations one can proceed as 
before to analyze the solutions. Depending on the precise parameters and the boundary 
value fo one finds again in the case of a semi-infinite chain regular or diverging solutions. 
The regular solutions can be periodic or quasi-periodic. In the case of a chain with finite 
size only a discrete subset of the above regular solutions survives. We are not aware of 
any tractable analytic method that determines these regular solutions for generic sets of 
parameters. 

4.2.1 Solitons and pinning of superconductivity 

It is known that the DNLS equation admits also another interesting kind of solutions: 
soliton and kink solutions. A detailed discussion of these solutions in the real domain 
and related references can be found in [B7J whose main points can be summarized briefly 
as follows. A priori one might expect that the DNLS equation does not admit such so- 
lutions. Soliton-like solutions are typically associated to integrable systems and DNLS 
is not integrable. It exhibits irregular chaotic behavior which in principle may prevent 
perfect localization. Nevertheless, it can be shown that non-integrability and discreteness 
appropriately combine to make such solutions possible. 

The solutions of interest have the following characteristics. They are solutions where 
the amplitude tp n is exponentially localized around a single site, say at n = 0. Following 



the nomenclature of [27] one can distinguish between two situations: 



(1) Bright solitons: in this case, \(p n \ > \(p n +x\ for n > and |<£> n | < |<^ n +i| f° r n < 
with limini^oo \<p n \ = 0. 

(2) Dark solitons: in this case \tp n \ < \ip n +i\ for n > and |y n +i| > Ifnl for n < with 
lim^i^oo \<f n \ > 0. It turns out that lim n _^ +00 ip n = — \im n ^_ OQ <p n , so these solutions 
are really kink solutions. 

In our context, where each site labeled by an index n, models a (1+1)- or (2+1)- 
dimensional layer of a superconducting material, such configurations would correspond in 
case (1) to situations where in a chain of layers the interlayer interactions work in such a 
way that energy and superconductivity are strongly localized around a central site. In case 
(2) the opposite happens. Energy and superconductivity are modulated in such a way that 
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they are almost uniform along the chain except around a central site where the condensate 
vanishes as it changes sign and superconductivity becomes very weak. 

In a continuum limit (see next subsection) the above configurations appear to recon- 
struct a junction of three materials with one dimension higher. In case (2) we recover a 
configuration that is very similar to the dark soliton of |T^] and reminds of an SNS junction 



of (2+1)- or (3+1)- dimensional superconductors. From this point of view the configura- 
tion of case (1) resembles a junction of two materials in the normal state separated by a 
thin superconducting layer in the middle. 

Since (p n are now real it is convenient to view the equations (|4.8|) as a two-dimensional 
map M. : M 2 — » M 2 by defining a new set of M 2 coordinates (x n ,y n ) = (<p n ,<p n -i). Then, 

M : ( Xn+1 = ~^ + ^ Xn ~ yn . (4.26) 

v Vn+1 X n 

The identification of the soliton-like solutions is closely related to the structure of the fixed 
points of this map. 

The fixed points, which by definition obey the relation x = y, are located at 



*o = 0, x ± = ±y-^±^. (4.27) 

The x± fixed points exist only when sgn(g + 2) = — sgn(s). Greene's residue p for the fixed 
point at the origin is [27| 



P=\(9 + 2). (4.28) 

Consequently, for \g\ < 2, we obtain < p < 1 which implies that the origin is a stable 
elliptic fixed point encircled by stable elliptic orbits. When in addition s < the fixed 
points x± are unstable hyperbolic fixed points. 

For \g\ > 2 we may distinguish between the following two cases: 

(i) g < —2, s > 0. In that case p < and the origin becomes an unstable hyperbolic 
point. The points x± are stable elliptic fixed points. 

(ii) g > 2, s < 0. The fixed point at the origin becomes unstable and through period- 
doubling bifurcation a new period-2 orbit appears located on the line x = —y. 

Before proceeding to explain the main idea underlying the existence of soliton-like 
solutions it will be useful to introduce some language which is common in the study of 
dynamical systems. 
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A set W is called an invariant manifold of a dynamical system if for any point x G W 
the dynamical evolution of x for any amount of time t continues to belong in W. Every 
fixed point p comes with its invariant manifolds. Such manifolds are called stable, and 
denoted as W s (p), if all points that belong on them approach asymptotically the fixed 
point p under dynamical evolution (namely p is an attractor on W s (p)). In contrast, an 
invariant manifold is called unstable, and denoted as W u (p), if all points that belong to 
it move asymptotically away from the fixed point p under dynamical evolution (in other 
words, p is a repellor on W u (p)). 

For generic non-integrable maps it is known that the stable and unstable invariant 
manifolds of hyperbolic fixed points cross each other. Points that reside on the intersection 
of stable and unstable invariant manifolds of the same fixed point are called homoclinic 
points. Accordingly, points that reside on the intersection of stable and unstable invariant 
manifolds of two different hyperbolic fixed points are called heteroclinic points. 

Having made this short introduction, we are now in position to describe what happens 
in our specific system provided by the map (|4.26|) . First consider the case (i) with g < —2, 
s > 0. The origin is an unstable hyperbolic point. Moving along an orbit on an unstable 
manifold W" of the origin and then crossing through a homoclinic point to a stable manifold 
W s gives rise to a soliton-like solution of the type (1) above. An explicit computation of 
the stable and unstable manifolds in the case of DNLS as well as specific examples can be 
found in [pTj] (see, for instance, Fig. 7 in f27| ). 

In case (ii) with g > 2, s < one obtains a similar soliton-like solution, but with the 
new feature that adjacent amplitudes <p n , <p n +i have alternating signs. Such solutions are 
known as staggered solitons E7\. The solitons of the previous paragraph are also known 
as unstaggered solitons. 

Finally, in the case of \g\ < 2, s < one can consider heteroclinic orbits connecting the 
two unstable hyperbolic fixed points x±. These are kink solutions of the type (2) above. 

4.3 Continuum limit and the Gross-Pitaevskii equation 

It is interesting to consider the continuum limit of the chain configuration described above. 



Th G S 

M °° ' Td ~* X ' ^ ~* ^ ' 9 IvP ' S ~* M2 ( 4 ' 29 ) 



In this limit 



with the new parameters x, G, and S kept finite. In this limit the recursion relations ( (4.5|) 
turn into the second order non-linear Schrodinger differential equation 

<P" + <P (g + ^\<p\*- 2 ) =0 (4.30) 
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where ' = -4-. 

ax 



For S = 4 we recover directly a well-known equation in the context of superfluids; the 

H which gives a coarse-grained description of superfluids 



Gross-Pitaevskii equation [48 



at long wavelengths. The GP equation is typically relevant for weakly-interacting Bose- 
Einstein condensates or strongly bound fermionic superfluids at low temperature. It is 
interesting that in our formalism, which has a radically different point of departure, the 
same description emerges naturally (as a candidate description of strongly coupled super- 
conductor physics) out of a framework designed specifically to deal with layered structures. 
It would be worth exploring further parallels that may exist between our formalism (and 
the generalizations of the GP equation that it suggests) and the known applications of the 
GP equation to superfluidity and superconductivity. 

For simplicity and concreteness let us continue to concentrate on the case of 5 = 4. 
Two well-known solutions of this equation are: 



(1) Bright solitons: for G < 0, S > 



<p(x) = ± 



4G 



S cosh (V— G x) 



(4.31) 



(2) Dark solitons: for G > 0, S < 



(p (x) — ±\ ——- tanh 
V ^ 




(4.32) 



For real (p we can find a more general class of solutions expressed in terms of the Jacobi 
elliptic function sn(u\m) 



, , . 2{G - VG 2 + SO) , ( 
ip(x) = ±i\j sn ( (x + xq) 



C and xq are integration constants. 

Reinstating the phase (|4.4| ) we obtain 



'G + VG 2 + SC 



a(x) = e ix0 <p(x) 



G - VG 2 + SC 



G + VG 2 + SC 



(4.33) 



(4.34) 



where 9 is the finite rescaled version of the angular inter-layer coupling d 



9 = Md 



(4.35) 
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Figure 5: Plot (a) depicts a generic sine-like periodic configuration for G > 0. The specific 
plot has G = S = 1. Plot (b) depicts a generic periodic configutation for G < 0, S > (in this 
particular case G = — 10, S = 1). By suitably tuning G,S one obtain the bright- soliton of eq. 
Q4.31 ), here depicted in plot (c) for G = — 1,S = 1. A dark soliton (or kink) solution can be 



obtained by tuning G, S so that the second argument in the Jacobi sine function becomes 1. An 
example of this case appears in plot (d) for G = —S = 1. 



At a boundary point Xf, the discrete equation (|4.25| ) becomes in the continuum limit 
( 4.29| ) a Dirichlet boundary condition 



(p(x b ) = . (4.36) 



By varying the parameters of the solution ( |4.33| ) one obtains qualitatively different 



behaviors. For concreteness, set C — 1, Xq = 0. For G > and generic S ^ one finds 
periodic sine-like solutions like the one depicted in plot (a) in Fig. |5|. For G < and S > 
the generic periodic solution, that appears as plot (6) in Fig. [5|, can be suitably tuned 
to obtain the bright-soliton solution ( |4.31| ) (see also plot (c) in Fig. |5|). For G > and 
S < one can find kink configurations by choosing the parameters G, S so that the second 
argument in the Jacobi sine function becomes 1. In that case, we recover the dark soliton 
solution ( |4.32p (an example of such a configuration appears in plot (d) of Fig. ||). 

By suitably truncating any of the solutions depicted in plots (a), (6), or (c) of Fig. 
[5] within an interval bounded by the location of two zeros of ip one obtains trivially a 
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finite-size Josephson junction chain in a continuum limit. For the plot (d) there is a single 
point where ip vanishes (the core of the kink solution), hence solutions of this type that 
are consistent with the Dirichlet boundary conditions (|4.36|) do not exist. 

4.4 Josephson current across a chain 

Following the discussion of subsection |3.2| we can write the total current across the chain 

(evaluated across two adjacent superconductors at positions n — 1, n) as a sum of two 
contributions 

Jtot = (Jn-l,») + ■ ( 4 - 37 ) 

J n -x,n is t ne current operator 

J n _i, n = i{O n ^O\ - 0\_ x O n ) (4.38) 

associated to the charge-transferring part of the [n — l,n) inter layer interaction in eq. 
( |4.1| ), and J^° s f^f son is the Josephson current associated to the backreaction of the system 
discussed in subsection |3]2|. The total current J tot is the total conserved current of the 
system, which, up to a potential multiplicative constant that we will keep implicit, equals 
the quantity J in eq. fl4.13|) . This is a site- independent quantity. In a chain with finite 
length the boundary conditions set J = and (J n -i,n) = — Jn-\ > n S ° n as 111 the two-site 
system of subsection |3.2|. More generally, however, 



J = (Jn-l,n) + J^tf S ° n . (4.39) 

Using the definition ( [4.38j ) we find that the current 

(J n -l,n) = - «„_!«*) = ^(e^^n-l^n _ e ~"Vn-l^n) ■ ( 4 - 40 ) 

In terms of the polar coordinates ( |4.9D we further obtain 

( J„_ lin > = -2r n _ x r n sin(^ + A9 n ) = -2 J^±^il (4.41) 



where A9 n and J were defined in eqs. ( 4.12|) and ( 4.13|) respectively. 



For a general solution of the recursion equations ( [4.3|) and $ ^ 0, 7r this is a link- 
dependent current. It is periodically or chaotically modulated in periodic or quasi-periodic 
solutions. In real solutions, where A9 n = mod ir, 

(J n -i,n) = ±2r„_ 1 r n sin . (4.42) 

For instance, in soliton-like solutions this current is very weak except around a central site. 
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Figure 6: An (un) conventional J J constructed as two semi-infinite AdS/CFT arrays of (2+1)- 
dimensional holographic superconductors linked to each other at a two-dimensional weak-link in- 
terface. 



In the special case where $ = 0, ir (i.e. when no external interlayer gauge field is 



applied), eq. ( |4.41|) gives the site- independent current 

(J„_ 1)n ) = ±2J . (4.43) 

This current vanishes for a chain with boundaries, but can be non-zero in chains without 
boundaries, e.g. in a circular chain with periodic boundary conditions. This is a simple 
example of how the topology of the network can affect the qualitative features of the 
configuration. 

4.5 Towards a typical Josephson junction 



In subsection |3.2| we discussed the similarities and differences between a two-site system 
and typical Josephson junctions. Here we discuss how the two-site system can be extended 
to look more like the typical Josephson junction. 

Assume we want to describe a junction of two superconductors in three spatial dimen- 
sions linked weakly across the third direction z at a two-dimensional interface. In previous 
subsections we described how to deconstruct (3 + l)-dimensional layered superconductors 
from an array of (2 + l)-dimensional holographic superconductors using linear AdS/CFT 
arrays. To construct an (un) conventional Josephson junction of two superconductors of 
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this type a possible strategy is described in Fig. ||. A layered superconductor on the left 
(right) is deconstructed as an array of cites linked through interactions of the form 

h L{R) (o^O L ™ + 0**»0«g) . (4.44) 

n 

With a real coupling hu^ no external transverse gauge field is applied along the z direction. 
Across the two-dimensional interface the right-most black site of the left chain can be linked 
to the left-most red site of the right chain through a link of a different type depending on 
the specific nature of the left and right sites. For s-wave holographic superconductors both 
on the left and the right a simple example of a double-trace weak link is 

hi ink {O l O R] + L] R ) . (4.45) 



Then one can solve the analog of the equations ( |4.3| ) and determine the Josephson current 
as was described in the previous subsection. SNS-type solutions of a uniform array with 
h>L — h,R — hu n k were described in subsections |4.1| , [4.2| , |4.3| (in a discrete or continuum 



limit). In general, the asymptotic difference of the phase of the condensates, Ai9 = '^l — '&r, 
is a dialed quantity in these systems. For conventional SNS or SIS-type JJs we anticipate 
the presence of a Josephson current that follows the sine law relation J max sin A$. We hope 
to return to a detailed survey of such systems in future work. 



5. Outline of future directions 

We have proposed a novel holographic way to model the physics of Josephson junctions 
using networks of (super) gravity theories on asymptotically AdS spacetimes coupled via 
mixed boundary conditions. One of the advantages of this approach, compared to previous 
holographic approaches, is the versatility by which it can incorporate many different types 
of Josephson junctions and networks with limitless possibilities in their architecture. For 
conventional SNS or SIS-type superconductors we presented a simple two-site model that 
exhibits some of the standard features of Josephson junction, e.g. the sine relation between 
the Josephson current and condensate phase difference. We explained in what sense this 
system is different from the typical Josephson junctions and how one can use AdS/CFT 
arrays to describe the more typical systems. We have also seen how a simple network on 
a chain produces complex dynamics with a variety of interesting features. 

Our preliminary analysis opens the possibility for a diverse set of calculations and 
extensions. Some of the most prominent ones are the following. 
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Figure 7: The architecture of the Y-Josephson junction network. 



(a) Finite temperature and a more complete analysis of phenomenological implications 

We have so far considered simple examples of holographic JJs at vanishing temperature 
and charge densities. It is of obvious interest to extend the setup to finite temperature 
using hairy black holes in designer multi-gravity (see e.g. for related work) and to 
explore possible phase transitions as we vary the temperature and/or charge densities. 

Extending the list of examples it is desirable to consider the explicit properties of 
other holographic JJs (or JJNs) built from different types of holographic superconductors 
(e.g. s-wave, or p-wave). For example, it will be interesting to define and study interlayer 
transport coefficients in such models. The ultimate goal is to explore the extent to which 
these constructions reproduce known phenomenological features of JJNs or layered super- 
conductor physics. For example, it would be interesting to reproduce previously observed 
non-sinusoidal current-phase relations in unconventional JJs. 



(b) Other network architectures and complex behavior 

It has been pointed by many authors (see e.g. ||26|| ) that the architecture of a JJN can have 
important implications on the physical properties of the system. It is interesting to explore 
other configurations and examine how they affect the collective and local properties of the 
sites. Networks with double-trace or higher multi-trace interactions can be constructed. 
An example that has been studied previously in the condensed matter literature is the 
Y-Josephson junction (see for instance ]5I| ). The architecture of a Y-Josephson junction 
network appears in Fig. [7]. 
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It may also be interesting to explore the existence of vortex solutions in two- or three- 
dimensional AdS / CFT lattices. This becomes even more interesting in view of the observed 
connection to the GP equation in the continuum limit. A recent discussion of vortex 
solutions to the GP equation (and a related AdS/CFT application from a different point 
of view) can be found in [J5TJ . 

(c) Continuous limits and deconstruction 



In subsection [O] we considered a continuum limit of a one- dimensional holographic JJN. 
In this limit the number of sites is scaled to infinity with an appropriate scaling of the other 
parameters of the system to zero. For a specific set of parameters we recover in this limit 
the Gross-Pitaevskii equation. It would be interesting to explore further relations between 
the generalizations of this equation suggested by our formalism and known applications of 
the GP methodology in superfluidity and superconductivity. 

In addition, it would be interesting to explore similar continuous limits of other JJNs 
with more complicated topology and different ingredients. Such limits may be used to 
simplify some aspects of the analysis of the network or in order to attempt a novel decon- 
struction of one or more extra spacetime dimensions (more comments on this aspect can 
be found in M. 
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Appendices 

A. On the linear stability of DNLS 

In this appendix we discuss in more detail the linear stability analysis of the DNLS equation 

(B 

g(p n + v3 n _i + (f n+ i + -(p n \(p n \ s ~ 2 = . (A.l) 
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Introducing a small perturbation u n around a solution ip n 



(o) 



(p n = v?i 0) + u n 



(A.2) 



we obtain (at first order) the equation 

«Wi + + (5 + f ki°Y- 2 ) «» + ! ^^ 0)2 |^° ) | 5 - 4 < = . (A.3) 

Next we decompose ipffl and u n into their real and imaginary parts (x n , y n here should not 
be confused with the corresponding variables in eqs. ( [4.14a ), ( |4.14b|) in the main text) 



n ^ X n ~\- iY n , U n X n -\~ i]Jr> 



<P 



For these variables we obtain the following two coupled sets of equations 
x n+ i + z n _i + (g + ^{Xl + Yl 
+£ (5-2) + ^ ((X 2 _ Y 2 )Xn + 2Xn Y n y n ) = 



2 X r 



(A.4) 



(A.5a) 



4 



Introducing the notation 



J/n+l + yn-1 + + + ) ^ ] 2/n 

+£ (5-2) ^ + (_ (X 2 _ Y 2 )yn + 2X n F n X n ) = 



5(5-2) 



— % 



we can rewrite the set of equations ( |A.5a ), ( |A.5b|) as a matrix equation 

/ x n+1 \ (Ml -1 N n \ / x n \ ( x n \ 

10 
N n M^ -1 
\ 1 J 

The characteristic polynomial for the eigenvalues A of the matrix J„ is 
A 4 - (M* + (1 + A 2 )A + A 2 [M x n Ml + 2 - X 2 ) + 1 = . 



Vn+l 
\ Vn ) 



■En— 1 



•En—1 
tin 



(A.5b) 

(A.6) 

(A.7) 
(A.8) 



(A.9) 



(A.10) 
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By further setting 



M» = -5 - tX + Y, 



we can recast ( |A.10|) into the more convenient form 

(A 2 -/inA + l) 2 = A 2 ^(X^ + r n 2 ) 2 
which yields four roots labeled by two Z 2 indices £1,62 = ± 



Hn + e^niXi + Y%) + e 2 \ (//„ + eny n (X^ + Y*)y - 4 



(All) 



(A.12) 



(A.13) 



The discriminant 

A = + e,v n [Xl + Y*)) 2 -4>=(g + l{8- 8,(5 - 2))(X 2 + F n 2 )^ * - 4 (A.14) 

can be either positive or negative. When A < 0, the eigenvalues A £lj£2 lie on the unit circle 
and the solution is linearly stable. Since 



A = A , A. 



(A.15) 



with 



A ± = ~g±2 + -(5-e 1 (6-2))(X 2 n + YZ 



A < requires (by definition A + > A_) 



A . > 0, A_ < . 



(A.16) 



(A.17) 



For s > the second term on the rhs of eq. ( A.16 ) is always positive (recall that 5 > 2 
for A < |). Requiring ( |A.17| ) gives 



s<2 , _^ <(M)¥< §|^. 

s s(o — 1) 

The lower bound on (X 2 + Y^f~^ is trivial when in addition —2 < g. 
For s < the analog of ( |A.18|) is 



2<~a 2{2 -~ 9) ,(X* I Y")^ c 2{ ~ 9 + 2) 



(A.18) 



(A.19) 
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